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O , Abstract 

Supported lipid membranes are useful and important model systems for studying cell 
> , 
^O i membrane properties and membrane mediated processes. One attractive application of 

supported membranes is the design of phantom cells exhibiting well defined adhesive prop- 

r — I , erties and receptor densities. Adhesion of membranes may be achieved by specific and 



non-specific interactions, and typically requires the clustering of many adhesion bonds into 
"adhesion domains". One potential mediator of the early stages of the aggregation pro- 
cess is the Casimir-type forces between adhesion sites induced by the membrane thermal 
fluctuations. In this review, I will present a theoretical analysis of fluctuation induced ag- 
gregation of adhesion sites in supported membranes. I will first discuss the influence of a 
single attachment point on the spectrum of membrane thermal fluctuations, from which the 
free energy cost of the attachment point will be deduced. I will then analyze the problem 
of a supported membrane with two adhesion points. Using scaling arguments and Monte 
Carlo simulations, I will demonstrate that two adhesion points attract each other via an 
infinitely long range effective potential that grows logarithmically with the pair distance. 
Finally, I will discuss the many-body nature of the fluctuation induced interactions. I will 
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show that while these interactions alone are not sufficient to allow the formation of ag- 
gregation clusters, they greatly reduce the strength of the residual interactions required 
to facilitate cluster formation. Specifically, for adhesion molecules interacting via a short 
range attractive potential, the strength of the direct interactions required for aggregation 
is reduced by about a factor of two to below the thermal energy ksT. 



1 Introduction 

Fatty acids and other lipids are essential to every living organism. Because of their amphiphilic 
nature, they spontaneously self-assemble into bilayer membranes that define the limits of cells and 
serve as permeability barrier to prevent proteins, ions and metabolites from leaking out of the cell 
and unwanted toxins leaking in [1] . In euokaryotic cells, membranes also surround the organelles 
allowing for organization of biological processes through compartmentalization. In addition, 
biological membranes host numerous proteins that are crucial for the mechanical stability of 
the cell, and which carry out a variety of functions such as energy and signal transduction, 
communication, and cellular homeostasis |2]. 

An important aspect of biological membranes is that they are typically not free but rather 
confined by other surrounding membranes, adhere to other membranes, and attach to elastic 
networks like the cytoskeleton and the extracellular matrix. Several model systems with reduced 
compositional complexity have been designed to mimic biological membranes. These biomimetic 
systems include phospholipid bilayers deposited onto solid substrates (solid-supported mem- 
branes) |3|, or on ultra-thin polymer supports (polymer-supported membranes) |1]. Placing a 
membrane on a flat substrate allows for the application of several different surface sensitive tech- 
niques, including atomic force microscopy, x-ray and neutron diffraction, ellipsometry, nuclear 
magnetic resonance, and others [5j. With the aid of biochemical tools and generic engineer- 
ing, supported membranes can be functionalized with various membrane-associated proteins [6]. 
Synthetic supported membranes with reconstituted proteins are increasingly used as controlled 
idealized models for studying key properties of cellular membranes [7] . They provide a natural 
environment for the immobilization of proteins under nondenaturating conditions and in well- 
defined orientations [8]. Another attractive application of supported membranes is the design 
of phantom cells exhibiting well defined adhesive properties and receptor densities [H]. Using 
advanced imaging techniques, detailed information can be obtained about the structure of the 
adhesion zone between the receptor-functionalized supported membrane and ligand- containing 
vesicles that can bind to the supported membrane [TUIITT] . These studies provide insight into the 
dynamics of adhesion processes and the molecular interactions involved in cell adhesion [T2|[T3]. 
Understanding these interactions is crucial for the development of drug delivery systems that 



depend on efficient adhesion between a liposome and the plasma membrane of the target cell. 

Adhesion between two membranes or between a membrane and another surface can, in princi- 
ple, be facilitated by non-specific attractive interactions (e.g.. Coulomb and van der Waals inter- 
actions) [HHT7] . Cell adhesion, however, is usually caused by highly specific receptor molecules 
located at the outside of the plasma membrane of the cell, that can bind to specific ligands 
on the opposite surface [IHIII^. Typically, the area density of the receptor molecules located 
at the outside of the plasma membrane is rather low which does not lead to efficient adhesion. 
However, when facing a surface with enough ligands, the receptors may cluster into highly con- 
centrated adhesion domains to establish much stronger binding [20l[2T]. Formation of adhesion 
clusters occurs in many biological processes [22], including the binding of white blood cells to 
pathogens [23j, cadherin-mediated adhesion of neighboring cells ^24j, and focal adhesion of cells 
to the extracellular matrix |25] . Many biophysical aspects of specific adhesion processes, ranging 
from the cooperativity in adhesion cluster formation to the influence of stochastic processes such 
as the ligand-receptor reaction kinetics, have been and continue to be studied theoretically using 
various models [261437] . 

Adhesion induced domain formation requires some attractive intermolecular interactions be- 
tween the receptor-ligand pairs. These interactions include both direct and membrane-mediated 
contributions. The former are typically described by pairwise potentials which are infinitely 
repulsive at very small molecular separations and attractive at somewhat longer (but still finite) 
distances |38]. Their effect can, therefore, be studied in the framework of the thoroughly re- 
searched lattice-gas model [32]. In contrast, much less is known about the membrane-mediated 
mechanism, which has been proposed by Braun et al. to explain to formation of gap junctional 
plaque at cell-cell interfaces [IQ], and whose origin can be understood as follows: Consider two 
adhesion bonds between two membranes or between a membrane and a surface (Fig. [Tj^A)). 
The adhesion points restrict the thermal height fluctuations of the membrane in their vicinity. 
This entropy loss can be minimized if the two adhesion bonds are brought to the same place 
(Fig. [U^B)), in which case the membrane becomes pinned at only one place rather than two. The 
membrane fluctuations, thus, induce an attractive potential of mean force between the adhesion 
bonds. This effect is often named after Casimir who predicted the existence of an attractive 



force between two conducting plates, due to quantum fluctuations of the electromagnetic field 
in the intervening space [H]. Later, Fisher and de Gennes generalized this concept to classical 
interactions induced by thermal fluctuations in soft matter systems |12]. For bilayer membranes, 
there is a great body of theoretical work on the Casimir effect between trans-membrane proteins 
(see review in [^, and refs. therein). Just like adhesion bonds, membrane inclusions represent 
a "constraint" on the shape of the membrane and, therefore, one should expects that they also 
interact with each other through Casimir-like interactions. In addition to the fluctuation-induced 
forces, the inclusions also experience other membrane-mediated interactions which arise from the 
membrane curvature elasticity and from the packing of the lipids near the inclusions' surfaces 
(see review in [Hj, and refs. therein). These other types of membrane-mediated interactions are 
also expected to exist between membrane adhesion bonds. 





Figure 1: (A) Schematic of a membrane attached by two distant adhesion bonds to an underlying 
surface. There is an entropy penalty associated with each adhesion bonds due to the restrictions 
imposed on the membrane thermal fluctuations in their vicinity. (B) The entropy cost can 
be minimized by bringing the adhesion bonds close to each other, in which case the thermal 
fluctuations become limited at only one location. The increase in the entropy in (B) compared 
to (A) is the origin of the attractive fluctuation-induced interactions between the adhesion bonds. 

The fundamental difficulty in attempting to provide a statistical-mechanical analysis of the 
aggregation behavior of the adhesion bonds is the need to integrate out the membrane degrees 
of freedom and write down the potential of mean force as a function of the coordinates of the 
adhesion sites 0(ri,r2,r3, . . . ,'r/v). This is a non-trivial problem since the membrane-mediated 
potential 0(^*1,^2, rs, . . . , rjy) is a many-body potential which cannot be expressed as the sum of 
two body terms. The many-body nature of 0(ri, "^2, "^3, • • • , "Hv) is best illustrated by the following 
example: Consider the configuration shown in Fig. [2]^ A) with two adhesion bonds at located at 
r*! and r2 and, in comparison, the one shown in Fig. [2|^B) with a single bond at n and a cluster 
of three bonds around r*2- Clearly, the spectrum of membrane thermal fiuctuations in both cases 
is quite the same and, therefore, the adhesion bond located at fi is attracted to the three-point 
cluster in[2](B) by the same force to which it is attracted to the single adhesion point in[2|^A). If 



0(ri, r2, rs, . . . , r/v) was the sum of pair interactions, the force in Fig. [2]^B) would be three times 
larger than the force in [2]( A) . 





7^2 n r2 ri 

Figure 2: (A) Schematic of a supported membrane with two adhesion located at fi and r2- 
(B) Similar to (A), but with a three-bond cluster instead of a single adhesion bond in r2- The 
adhesion bond in n is equally attracted (by a Casimir-like force) to the adhesion bond located 
in r2 in (A) and to the cluster of three adhesion bonds shown in (B). 



2 Lattice-gas model for adhesion bonds in supported mem- 
branes 

What is the difference between the aggregation of adhesion bonds in supported membranes and 
the traditional process of gas to liquid condensation? Condensation phase transitions are usually 
associated with a competition between the mixing entropy S which is higher in the dilute gas 
phase, and the interaction energy U which is lower in the condensed liquid state. The equilibrium 
phase corresponds to the minimum of the free energy F = U — TS, where T is the temperature 
of the system. At high T, the free energy F is "entropy-dominated" and equilibrium is attained 
in the gas phase. Conversely, at low T, the free energy is "energy-dominated" and, therefore, the 
condensed phase becomes thermodynamically more favorable. The liquid-gas phase transition 
can be analyzed in the framework of an Ising-like model of identical particles that populate a 
lattice. Excluded volume interactions between the particles are represented by the fact each 
lattice site can be occupied by no more than one particle. When two particles occupy nearest- 
neighbor sites, they interact in a pairwise fashion with an attractive energy — e. Denoting the 
occupancy of a lattice site by Sj, with Sj = for an empty site and Sj = 1 for an occupied site. 



the Hamiltonian of the lattice-gas model is given by 

y-LG = -e'^SiSj, (1) 

where the sum runs over all the pairs of lattice nearest neighbor sites. The phase diagram of the 
lattice-gas model is well known. There exists a critical value Oc such that if the interaction energy 
e < acksT, the particles will be distributed uniformly within the lattice. Above this critical value, 
e > acksT, a uniform distribution of the particles is observed only at low concentrations of 
particles ("gas phase"), but upon increasing the concentration of particles, the system undergoes 
a first order phase transition and a second coexisting phase appears with a considerably larger 
concentration ("condensed phase"). 

As discussed in the previous section, the aggregation process of adhesion domains involves an 
additional attractive potential of mean force resulting from the membrane thermal fluctuations. 
A lattice model where each lattice particle represents an adhesion bond in a supported membrane 
can, therefore, be introduced by supplementing Eq. ([1]) with an energy term corresponding to the 
fluctuation-induced interactions. Since the functional form of this many-body potential of is yet 
unknown, we would, at this moment, introduce it via a general potential function (p that depends 
on the coordinates of the lattice particles: 

n = -eJ2s,s,+<l)i{s..}). (2) 

Our first task must be to derive an expression for ({sj}). Once this is accomplished, one can 
attempt to analyze the statistical mechanical properties of the model and address the question 
appearing at the beginning of section |2l One particular issue that we would like to address is 
whether the fluctuation induced attractive potential (which is of entropic origin) can win the 
competition against the repulsive force originating from the mixing entropy? In other words, 
can adhesion clusters form for purely entropic grounds, i.e. for e = in Eq.([2])? Gas to liquid 
condensation transitions are generally believed to involve energy vs. entropy competition [15], 
but purely entropic phase transitions from a fluid (disordered) phase into a crystalline (ordered) 
phase are known to exist. Hard sphere systems, for instance, undergo a first order phase transition 



from a low density fluid phase into a high density solid phase [16]. This transition results from 
the competition between two entropies - the configurational mixing entropy which is higher in 
the disordered phase, and the entropy associated with the free volume available for each sphere, 
which is higher in the ordered crystal. 



3 Statistical mechanics of a membrane with one adhesion 
point 

We start our analysis by considering the system shown schematically in Fig. |3l consisting of a 
membrane with bending rigidity k that fluctuates above a flat impenetrable surface jl^. Let 
h (f) > be the height of the membrane with respect to the surface, and assume that the 
membrane is pinned to the surface at one fixed point located at ro {h (fo) = 0). The elastic 
curvature energy of the membrane is given by the Helfrich effective Hamiltonian 
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(3) 



where $ represents the hard wall constraint due to the surface ($ = 1 for h > 0, and $ = +oo 
ioT h < 0), 6 is the Dirac delta- function, and the integration is taken over the cross sectional 
(projected) area of the membranes of size L^. 




Figure 3: (A) Schematic picture of the model system consisting of a membrane that fluctuates 
above a flat impenetrable surface to which it is pinned at a single point. 



One can calculate the partition function Z corresponding to Hamiltonian (^, by considering 
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Figure 4: (A) A freely fluctuating membrane can be always translated vertically such that the 
point at which its height function h [r) attains its global minimum is on a flat surface and the rest 
of the membrane is above the surface. (B) A freely fluctuating membrane can be also translated 
horizontally. All the membrane conflgurations generated in this way are similar to each other, 
and the one for which the point of absolute minimum is at fo (represented by the solid line) is 
identical to the pinned membrane conflguration shown in Fig. |3l 



the Helfrich effective Hamiltonian of a freely fluctuating membrane 
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In this case, the associated partition function Zfrcc is readily calculated by introducing the Fourier 
transformation of h{r): hq = {1/L'^) J h{f)exp{iq ■ f), which decompose Hamiltonian (j4]) into 
the sum of independent harmonic oscillators 
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where / is a microscopic length scale of the order of the bilayer thickness. Hamiltonian (|3]) 
which also includes the functions $ and 6 cannot be diagonalized in the same manner. However, 
one can relate the partition function Z of Hamiltonian ([3]) with the partition Zfree of the free 
membrane Hamiltonian (jlj), by using the following simple argument. The energy of a freely 
fluctuating membrane is invariant with respect to rigid-body transformations such as a vertical 
translation {h (r) — )■ h{f) — Hq) of the membrane's center of mass. Therefore, one can draw 
a flat surface and translate the free membrane such that the global minimum of its height 
function coincides with the surface (see Fig. 111(A)). The vertically translated free membrane 
looks very similar to the pinned membrane shown in Fig. [3]^ A). The only difference between 
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them is that the former can also ghde over the surface (see Fig. Hl^B)), while the latter is pinned 
at a fixed position on the surface. This suggests that the pinning point effectively eliminates the 
membrane horizontal translational degree of freedom. In a statistical mechanical language, the 
configurational phase space of the pinned membranes is smaller than, yet similar to, the phase 
space of a free membrane. Each subspace of identical free membrane configurations, like the 
ones shown in Fig. IH^B), includes one pinned membrane configuration - the configuration where 
the minimum of h (r) is at the pinning point tq (or, more precisely, within a microscopic area 
of size P around the pinning point, where / is the spatial resolution of the continuum model). 
This pinned membrane configuration occupies a fraction {l/L)"^ of the corresponding larger free 
membrane configurational subspace, which implies that the partition functions of the two systems 
are related by Z = {1/ LYZ^^^c- The free energy is obtained from 

F = -kBT\n{Z) = -ksT In (Z^ee) + ^ksT In (j^ (6) 

The first term on the right hand side is the free energy of the free membrane whose elastic energy 
is given by Helfrich Hamiltonian ([3]). The second term, 

-^attachment,! = 2kBT lu ( — 1 , (7) 

is the free energy cost of attaching the membrane to the surface at one point. 

Following the above argument leads to a very interesting conclusion. Because of the similarity 
mapping that exists between the configurational phase spaces of the two problems, the statistical 
properties of the pinned and the free membrane must be identical to each other. This surprising 
result can be demonstrated by using an implicit-solvent coarse-grained (ISCG) bilayer model 
which enables molecular simulations of mesoscopically large bilayer membranes over relatively 
large time-scales [IHHSI]- Towards this end, we ran two independent Monte Carlo (MC) simula- 
tions - one of a free membrane (without a surface) and one of a membrane supported by a flat 
impenetrable surface. A snapshot from the supported membrane simulations is shown in Fig. [51 
Each lipid molecule is represented in the model by a short string of three spherical beads, where 
one of the beads (depicted as a dark gray sphere in Fig. [5]) represents the hydrophilic head group 
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Figure 5: Equilibrium configuration of a membrane consisting of 2000 lipids. Each lipid is 
represented by a trimer of one "hydrophilic" bead (dark gray sphere) and two "hydrophobic" 
beads (light gray spheres). The membrane is fluctuating above a plane surface (frame indicated 
by a thick black line), while one of the hydrophilic beads (the black sphere appearing at the front 
of the figure and indicated by an arrow) is held on the surface at a fixed position. 

and two beads (light gray spheres in Fig. [5]) represent the hydrophobic tail of the lipid. In the 
supported membrane simulations, the head bead of one of the lipids (appearing in the corner at 
the front of the figure and indicated by an arrow) was fixed to a flat surface which the lipids 
were not allowed to cross. We measured the Fourier spectrum of the membrane height function. 
For the free membrane, the application of the equipartition theorem to the Fourier-space rep- 
resentation of the Helfrich Hamiltonian (jS]) yields the following relationship between the mean 
squared amplitude of the Fourier modes ( "spectral intensity" ) and the wave- vector q: 

M - '-0- (8) 

Fig. E] depicts the results of our MC simulations for the spectral intensity vs. the wavenumber n = 
qL /{2tx). The figure show that, in agreement with our predictions: (i) the free (open circles) and 
pinned (solid circles) membranes exhibit the same statistics of thermal height fluctuations, and 
(ii) the spectral intensities of both membranes follow the tT'^ power-law dependence anticipated 
by Eq.([8D (dashed line). 

One can reverse the argument and derive Eq. ([7]) starting from the assumption that the spec- 
tral intensity of the supported membrane is identical to that of a free membrane and, therefore, 
can be described by Eq.(|H]). The derivation proceeds as follows: First, from Eq.(|S]), one can 

11 



CM 



CM 



A 



V 



1 


•n 


1 

^®. 

•>> 

1 


1 

"S 
1 


1 1 


1 1 1 1 1 1 


: 


0.1 
0.01 




o free membrane 

• supported membrane 

-4 

~ n 


- 


1 1 


1 1 


- 



10 



n 



Figure 6: The mean square amplitude of the thermal height fluctuations as a function of the 
wavenumber n. The results from the supported membrane simulations are shown by small 
solid circles. These results are essentially identical to those obtained from simulations of a free 
membrane which are represented by larger open circles. The dashed line indicates the asymptotic 
(|/iqp) ~ n~^ power law for small n. (adapted from [17] ) 
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Figure 7: The fact that the statistics of thermal height fluctuations is not affected by the single 
pinning point implies that the typical height of the fluctuations scales linearly with the distance 
from the pinning site. 

straightforwardly show that the typical height at which the membrane undulates above the 
surface at a distance r away from the pinning point scales linearly with r [171 152] (see Fig. [7]): 



u{r) = {h (r)) ~ r 



ksT 



K 



(9) 
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There is a repulsive force acting between the fluctuating membrane and the underlying sur- 
face, caused by their mutual steric hindrance. Helfrich [53] showed that the associated repul- 
sive interaction free energy density (per unit area) has the following scaling behavior V{r) ~ 
{ksTy / Ku{r)'^ which, together with Eq.dH]), yields 

V{r) ~ M:. (10) 

By integrating this energy density over the projected area of the membrane, one derives Eq.([7]) 
up to a numerical prefactor 

^attachment,! = / V{r)£rr^ f 2nr^dr = CksT In (j\ . (11) 

To set C = 2, as in Eq.(I7]), one needs to replace the scaling relation Eq.f jTOj) with the equality 

nr) - i^. (12) 



4 Fluctuation induced attraction between two adhesion 
points 



As noted by Helfrich [53] , the free energy density Eq. ( 1121) due to the steric hindrance between the 
two surfaces (i.e., the fluctuating membrane and the underlying supporting surface) is directly 
related to the rate of collisions between them. In other words, the probability density that the 
membrane hits the supporting surface at a distance r from the pinning point exhibit the same 
dependence on r as V{r): 

p[h{r) = 0]^^. (13) 

This relationship provides the information needed for calculating the fluctuation induced attrac- 
tive potential between two adhesion points. This is done by regarding the point of collision 
between the membrane and the surface as a second pinning point which can diffuse across the 
surface. In this context, the probability density p[h{r) = 0] is identifled with the pair correlation 
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Figure 8: The pair correlation function, g{r), of a non-stressed membrane vs. the pair distance 
r. The slope of the dashed straight line is —2. (adapted from [51]) 

function between the adhesion points which, therefore, also follows the scaling form 



9[r) 



(14) 



By definition, the pair potential of mean force is given by 



(f) = -keT In [g (f)] = 2kBT In (r) 



(15) 



which is an infinitely long range attractive potential that does not depend of the bending rigidity 
of the membrane, k. 

The validity of Eq. flT^ can be tested by using MC simulations of the ISCG model shown in 
Fig. [5] with two lipid heads attached to surface - one fixed at the origin and the other allowed to 
diffuse on the fiat surface. The pair correlation function is then directly computed by sampling 
the position of the mobile adhesion point. Our results [M], which are shown in Fig. [HI agree 
very well with Eq. (fT^ . The slope of the straight line on the log-log plot is equal to —2. The 
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Figure 9: (A) The typical height of the fluctuations of a supported membrane experiencing a 
harmonic confining surface potential, grows linearly close to the pinning point and saturates 
at large distances. (B) The pair correlation function, 5'(r), of a such a membrane vs. the pair 
distance r. The slopes of the solid and dashed straight lines are —2 and 0, respectively. ((B) is 
adapted from [51]) 

deviations from the power law behavior g{r) ~ 1/r^ at small values of r {r/L < 0.05) are related 
to the breakdown of the continuum description of the Helfrich Hamiltonian at small spatial scales. 
At small separations, the molecular nature of the lipids becomes important and the radial pair 
distribution function is dominated by the depletion shells around the lipids. 

What if, in addition to the excluded volume repulsion, the membrane and the surface also 
interact via an attractive potential of somewhat longer range? Let us consider, for instance, the 
case when a harmonic confining potential is added to the Helfrich Hamiltonian: 
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For the harmonically confined membrane, one can define the length scale ^^ ~ (K/7)^/^ which 
marks the transition between two scaling regimes. For r ^ C,j, the thermal fiuctuations are 
governed by the bending rigidity term in Hamiltonian (ITB]) . while for r ^ ^^ the harmonic 
confinement term becomes dominant. The latter term is a local one, which implies that the 
infiuence of the adhesion point becomes screened at large distances. In the case of a single 
adhesion point, the height of the fiuctuations is now given by (compare with Eq.([9])) 



, rJ^ forr<e7 

(Mr))- V^ ^^ , (17) 

U'-^ forr»e. 



as illustrated schematically in Fig. [HI^A). The correlation function of a pair of adhesion points is 
given by (compare with Eq. flT^ ) 

^-2 £qj^ r <^ ^^ 



9{r) 

57 



r° for r ^ ^^ 



The results of MC simulations of an ISCG molecular model of a harmonically confined membrane 
verify this crossover between the two scaling regimes of g{r) (see Fig. EJ^B)). 

The energy of a membrane subjected to lateral surface tension a > is given by the following 
Hamiltonian 



^-Ili(^'"y-^iH' 



d^f. (19) 



Scaling arguments [53] show that, in this case, the pair correlation function exhibits behavior 
intermediate between Eqs.(fT4l) and ( fTSl) : 

9 (r) -{ ' , (20) 

r ^ for r ^ No- 
where the crossover length ^o- ~ (k/ct)^/^. This scaling form is also confirmed by MC simulations 
(see Fig. [TO]). 
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Figure 10: The pair correlation function, g {"r), of a supported membrane under tension vs. the 
pair distance r. The slopes of the solid and dashed straight lines are —2 and —1, respectively, 
(adapted from [5^ ) 

5 The strength of the fluctuation-induced attraction 

One of the main questions we aim to explore is whether the fluctuation-induced interactions are 
sufficiently strong to allow the formation of adhesion clusters. In the case of two adhesion points 
the answer is no. Despite of the attractive force between the adhesion points, they remains 
unbound and their mean pair separation grows linearly with the system size L. More generally, 
if the pair correlation function decays algebraically at large distance, g {r) ~ r~'^, the mean pair 
separation is given by 



(r) 



^^^2 



Ji r'^gir) dr 
Ji f^gi"^ dr 



L for c < 2 

L/ In L for c = 2 

L^-" for 2 < c < 3 

In L for c = 3 

/ for c > 3 



(21) 
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The physically relevant cases in Eq. (!2T]) are c = 2, 1, and which, respectively, correspond to 
pinned, pinned-stressed, and pinned-confined membranes. In all of these cases, (r) grows with 
the size of the system. 

Another quantity of interest is the mean number (C) of contacts between the surface and 
a membrane with one adhesion point. As discussed in section |H the probability density that 
membrane comes into contact with the surface at a distance r from the pinning point has the 
same scaling form as the pair correlation function g (r). Thus, 



{C) 





{L/Q^ fore 


g{r)rdr ~ < 


L/^„ for c 




ln(L//) for c 





1 • (22) 

2 



We can use this last result to generalize and recalculate the attachment free energy of one ad- 
hesion point, Eq.(I71). Our original derivation of Eq.(I7]) was based on the argument that the 
configuration phase space of a pinned membrane comprises a small subspace within the configu- 
ration phase space of a free membrane. More precisely, we argued that this subspace includes the 
free membrane configurations in which the global minimum of the height function occurs at the 
pinning point of the corresponding supported membrane. We further argued that the relative 
size of the subspace is (l/L)"^, which was based on the assumption that typically there is only 
one contact point with the surface and, therefore, this contact point must be the adhesion site. 
However, as we see from Eq.( l22l) . a typical configuration makes (C) contacts with the surface. 
Therefore, the partition functions of the two problems (free vs. pinned membranes) are actually 
related hy Z = [(C) (l/L) ] Zf^cc- The attachment free energy is given by 



attachment,! — "n-B-t ili I — 7y 



2kBT\n{^^/l) forc = 

kBT\n{L/l) + kBTln{^Jl) for c = 1 • (23) 

2kBT \n{L/l) - kBT In [ln(L//)] for c = 2 



Notice that for sufficiently large L, Fattachment,l(c = 2) > Fattachmcnt,l(c = 1) > -Fattachment,l(c = 0). 

Indeed, the attachment of a free membrane to a surface is likely to be more costly than the 
attachment of stressed and harmonically confined membranes that exhibit reduced fluctuations 
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and, thus, remain close to the surface anyway. 

6 The many-body problem 

Let us look back at Fig. [8] which shows the pair correlation function between two adhesion 
points. The figure demonstrates that the scaling form Eq. (1141) holds over almost the entire range 
of pair separations considered {I < r < L/\^). The deviations from the power law at small 
pair distances arising from the short range depletion forces between lipids have already been 
discussed in section HI What is quite surprising, though, is the pretty good agreement between 
the MC results and Eq.( TT4l) at large pair distances. Eq.l IT^ has been derived for two adhesion 
points in a very large membrane, neglecting boundary effects. In the simulations the conditions 
are different - the membrane has a finite size and periodic boundary conditions are employed 
to reduce the finite size effects. Thus, each adhesion point interacts not only with the other 
adhesion point but also with its infinite array of periodic images. Nevertheless, the existence of 
periodic images seems to have a very small impact on the results. This observation is particularly 
unexpected for r > L/2 corresponding to situations where one of the adhesion points is equally 
close to two images of the other adhesion point. The only possible way to explain this surprising 
observation is to assume that the periodic images of the adhesion points are largely screened. 
This assumption is consistent with the following physical picture: The membrane mediated 
interactions originate from the entropic cost due to the suppression of the membrane thermal 
undulations. Thus, the presence of each adhesion point is felt only in the region where it affects 
the fluctuations and cause their reduction, while outside of this region, the adhesion point is 
effectively screened. In this perspective, the idea that distant adhesion points are screened seems 
logical. The fluctuations vanish at each adhesion point, irrespective of the distribution of the 
others. Moreover, in the immediate vicinity of each adhesion point, one expects the amplitude 
of the fluctuations to depend only on the distance from that adhesion point. If the membrane is 
neither stressed nor experiencing a conflning surface potential, the amplitude of the fluctuations 
in this region grows linearly with the distance r from the adhesion point, as given by Eq.Q. We 
now wish to introduce a more general expression that holds over the entire area of the membrane 
and coincides with Eq.dH]) close to every adhesion point. Our suggestion is as follows [55J: In 
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each unit area of the membrane, the mean height of the membrane above the surface is given by 
(compare with Eq.(|9])) 

(/i(r))~rf,ni„y^' (24) 

where d^i-^ is the distance of the unit area from the nearest adhesion point. We also replace r 
with c/min in Eq. flT2|) for the attachment free energy density, which now reads 

V(r) = Itf. (25) 

The total attachment free energy of a given distribution of adhesion points is obtained by in- 
tegrating the attachment free energy density Eq.( l25l) over the entire membrane area. This 
calculation is done by constructing the Voronoi diagram of the distribution of adhesion points, 
integrating the free energy density with each cell (where in each cell the distance is measured 
from the adhesion point located in the cell, and a small region of microscopic size / around the 
point is excluded from the integral), and summing the contributions of the different cells: 

^attachment = XI / ~V ^^^ (26) 

In a lattice-gas model, the discrete analog of this expression applies 

I ~, — I (1 ^ '5j), (27) 



'rnm 



where the sum run over all the empty lattice sites (sj = 0) and /^ is the area per lattice site. 

As discussed in section El our main goal is to develop and use a lattice-gas model for the 
aggregation problem of adhesion points. In the model, each lattice point represents an adhesion 
point between the membrane and surface. The energy of a given configuration of lattice points 
is given by Eq.([2]), where the first term represents the short-range attraction between adhesion 
points and the second term is a many-body fluctuation-induced potential (j) ({sj}). Our journey 
to derive an expression for (j) ({sj}) started in section [3l and has finally reached the end. ({sj}) 
is a potential of mean force which, for a given distribution of adhesion points, is determined by 
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Figure 11: Schematic of a square membrane of linear size L with two adhesion points located at 
[x, y) = (±r/2, 0). The dashed line shows the boarder between the Voronoi cells of the adhesion 
points. 



tracing over all the relevant membrane configurations and calculating the free energy penalty 
associated with the reduced thermal fluctuations. Eq. (!27|) provides this expression by assigning 
a free energy cost with each empty lattice site that represents a fluctuating unit area of the 
supported membrane. Introducing Eq. (p7|) into Eq.(j2]), yields the energy function of our lattice- 
model of adhesion points 



H 



Yl ^^^^ + Yl 



knT ( I 
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6.1 The two-body problem revisited 



Let us see how one can re-derive Eq. fllSI) for the pair potential of mean force by calculating the 
attachment free energy Eq. fl2^ . Towards this end, consider the membrane shown schematically 
in Fig. [TT]with two adhesion points, each of which located a distance r/2 from the center of the 
membrane. The dashed line shows the boarder between the Voronoi cells of the adhesion points, 
where each cell extends over half of the area of the membrane. For the configuration shown in 
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Fig. [m the attachment free energy Eg. (1261) reads: 



L/2 

-rattachment,2 ~ ^ I '^V 
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Integrating over y yields, 
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(30) 



Assuming that I < r <^ L, the inverse tangent function in Eq. flSDl) can be approximated by the 
constant value of 7r/2 over most of the integration range. With this approximation, one gets 



i"attachment,2(r,L) ~ 2A;bT lu ( y j + 2A;BTln (y 



attachment,! 



(L) + 0(r). 



(31) 



The first term in eq.( l3T|) is the free energy cost of a single adhesion site [Eq.([7])], which is the 
expected value when the two adhesion points coincide (r ~ I) to form a single cluster. The 
second term, which represents the additional free energy cost associated with the separation of 
the adhesion points, is identified as the fluctuation induced pair potential, in agreement with 
Eq.([l5]). 

6.2 Mean field theory 

We now come back to the many-body problem and start with a mean field analysis of our lattice 
model Hamiltonian ( l28l) . Let us consider a lattice of Ns sites of which N < Ns sites are occupied 
by adhesion points. Let us further assume that the adhesion points form N^ < N adhesion 
clusters. The free energy of system includes three contributions: (i) the mixing entropy of the 
adhesion clusters, F^i^, (ii) the lattice-gas energy, Ej^q, of the direct interactions between the 
adhesion points [first term in Eq. (l28p ]. and (iii) the attachment free energy, Fn [second term in 
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Eq.( l28|) ]. The first free energy contribution is given by 



F ■ 
knT 



In. .. 



+ 2H11- P^' 



where c is the second virial coefficient. On average, each cluster consists of (N/Nc) adhesion 
points; and if we assume that it has a roughly circular shape than c ~ 4{N/Nc). Denoting the 
number densities of the adhesion points by p = N/N^, and of the clusters by p* = N^/Ng < p, 
the free energy of mixing per lattice site is given by 

^"^'^ p*[\n{p*)-l] + 2pp*. (33) 



NsknT 



The second contribution to the free energy is due to the direct interactions between the 
adhesion points. The ground state of the interaction energy Ei^q is achieved when a single 
circular adhesion domain with minimal surface is formed. If we set the ground state as the 
reference energy, the energy of an ensemble of clusters can be estimated as being proportional 
to the total length of the domain boundaries. For Nc circular clusters of size (N/Nc) we have 

^Tr\ Tr = ^ypp^ (34) 



NgkBT N, V N, 
where A, the associated dimensionless line tension, is proportional to the interaction energy e 

A = 2^/^Be, (35) 

and B is the mean number of nearest-neighbor vacant sites per occupied site on the boundary 
of a cluster (B — ;■ 1 for very large clusters). The sum of free energy contributions (133|) and (!M|) 
constitutes the total free energy density (per lattice site) of a 2D lattice gas of clusters: 

^"-"^ p* ln(p*) -p* + 2pp* + AVpp^. (36) 



NgkBT 

The third contribution to the attachment free energy can be estimated as follows. The clusters 
form Nc Voronoi cells, each of which has on average an area of A^^^ = {Ns/Nc)!"^. The attachment 
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free energy of each Voronoi cell is given by an equation similar to Eq.(I7]) for the attachment free 
energy of one adhesion point, but with A^i instead of the total membrane area L^. Thus 



^Ar = iV, 



^.Tln(| 



(37) 



and the attachment free energy density is given by 

1^ = -/^*1^(/'*)' (38) 

NsKbT 

which eliminates the first term in the lattice- gas free energy density [Eq. (!36|) ]. yielding 

^ ^'^'^ +-^^ = -p* + 2pp* + Xy^. (39) 



N^ksT N^ksT N^keT 

We consider a low density of adhesion sites p ^ 1 , which also implies a low number density of 
adhesion clusters since p* < p. By minimizing the free energy density we obtain the equilibrium 
value of the p* for the standard lattice-gas model [Eq. (15^ ] and for the adhesion points of a 
fluctuating supported membrane [Eq.f pQ]) ]. In both cases, the system undergoes a first order 
phase transition at Ai(p) from the gas phase (p* = p) to a condensed phase consisting of only a 
few clusters (p* ~ 0). Also, in both cases, the free energy reaches a maximum at intermediate 
densities (0 < p* < p). This free energy barrier for condensation disappears at the spinodal 
point A2(p) > Ai(p). For the lattice-gas problem we find 

A^G = l-2p-ln(p) 

A^^ = -4p-21n(p), (40) 

while for the adhesion points of fluctuating membranes we have 

Ai = l-2p 

As = 2-4p = 2Ai. (41) 

The results of Eqs. 04Up and 04 ip are summarized in Fig. [TWA) and [T^ B). respectively. The 
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Figure 12: The phase diagram of the adhesion sites calculated within the mean field approx- 
imation. (A) Eg. (140 p for the standard 2D lattice-gas model. (B) Eg. (1411) for adhesion points 
of fiuctuating membranes. Ai and A2 represent the first-order transition and spinodal lines, 
respectively, (adopted from |5"5j ) 

important points in the results are that: (i) Ai > 0, which means that the fluctuation induced 
interactions alone are not sufficient to induce aggregation of adhesion domains, but (ii) they 
greatly reduce the strength of the direct interactions reguired to facilitate cluster formation 
since Ai < A^*^' (and also A2 < Ag^). Below, we support these conclusions with MC simulations 
and show that for adhesion points of fluctuating membranes, the site-site cohesive energy e for 
the onset of aggregation falls below the thermal energy ksT. 

6.3 Monte Carlo simulation 

To further investigate the aggregation behavior in supported membranes, we performed MC 
simulations of both our lattice model of adhesion points and of the standard 2D lattice-gas 
model [55]. We simulated the system at two different densities p = N/Ng = 0.05 and p = 0.1, 
and for various values of e ranging from to 3 ksT. Snapshots taken from simulations for 
e = IksT and p = 0.1 are shown in Fig. [121 Fig. [T^ A) shows the initial configuration where the 
points are randomly distributed on the lattice. Figs. [TST B) and[T3l(C) show, respectively, typical 
eguilibrium configurations of the standard lattice-gas model and of our model of adhesion points. 
One clearly sees that for the same strength of the interaction energy e = Ifc^T, the standard 
lattice gas model remains in the gas phase, while the adhesion points (that, in addition to the 
direct interactions, also attract each other via the fiuctuation-induced mechanism) condense into 
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Figure 13: (A) Initial configurations of tfie simulations in which the points are randomly dis- 
tributed on the lattice. (B) Representative equilibrium configurations of the standard lattice gas 
model for p = 0.1 and e = Ifc^T. (C) Representative equilibrium configurations of our lattice 
model of adhesion points for the same values of p and e as in (B). 

a large cluster containing almost all the adhesion points. 

In order to determine the onset of the gas to liquid transition, we measured the average 
number of clusters in the system (where a cluster is defined as a set of neighboring occupied 
sites), and the mean value of the energy of direct interactions between sites, (-Elg) [first term 
in Eq. (!28|) ]. Our results are summarized in Fig. [TW A) (for p = 0.05) and [TW B) (for p = 0.1). 
For each p, we measured these quantities both for the standard lattice-gas model (open symbols 
and dash-dotted lines in Fig. [T^ and for the adhesion points model (solid symbols and solid 
lines in Fig. [T^ . The number of clusters is denoted by squares (values on the right y-axis of 
the figures), while (-Elg) is represented by circles (values on the left y-axis). The gas phase is 
characterized by a large number of small clusters, some of which may be of the size of a single 
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Figure 14: Left y-axis: The energy of direct interactions between sites, (-Elg); as a function of 
e, for (f) = 0.05 (A) and = 0.1 (B). Solid circles - results for our model for adhesion points. 
Open circles - results for the standard lattice-gas model. Right y-axis: The number of clusters 
as a function of e, for = 0.05 (A) and = 0.1 (B). Solid squares - results for our model for 
adhesion points. Open squares - results for the standard lattice-gas model, (adapted from [55J) 

site. Furthermore, since each occupied site has a relatively small number of neighboring occupied 
sites, the mean configurational energy (— -Elg) is relatively low. Conversely, when the sites form 
large clusters in the condensed phase, (— -Elg) is high, and the total number of clusters decreases 
(and in many cases, especially for large values of e, we simply observe only a single cluster in 
our system). Fig. UM exhibits an abrupt, clearly first-order, transition from a gas phase with a 
large number of clusters and small (— -Elg) to a condensed state with a small number of clusters 
and large (— E'lg)- The estimated values of e at the transition are (see vertical lines in Fig. UM : 
et — O.TfcfiT (p = 0.05) and e^ ^ O.GbksT (p = 0.1). In comparison (see also Fig. [H]), for 
the standard lattice-gas model, the transition values are roughly twice larger than these values: 



^LG 



lAbkeT (p = 0.05) and e^^ ~ IMbT (p = 0.1). 



Our computational results which show that the fluctuation mediated interactions reduce 
the strength of ej, are in a qualitative agreement with the mean field theory prediction. To 
make a quantitative comparison between the theory and the simulations, one needs to estimate 
the parameter B appearing in Eq. (l35p . Several reasons make such an estimation difficult and 
inaccurate: First, our non-standard mean field theory is based on the assumption that the 
clusters are circular and roughly have the same size, which is quite a crude approximation. 
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Second, tracing the precise location of et in Fig. [H] is largely inaccurate because of the finite size 
of the system that makes the transitions look like crossovers. To reduce the large uncertainties 
associated with the determination of et, one can look at the difference between the value of this 
quantity in our model of adhesion points and for the standard lattice-gas model. Using 

A^« - Ai = 2V^B (e,L« - et) , (42) 

for p = 0.1, we find B c^ 1, as indeed expected for large clusters. 

7 Conclusions 

In this review we presented a statistical thermodynamics analysis of the aggregation behavior of 
adhesion points between a fluctuating membrane and a supporting surface. Our analysis focused 
on the contribution of the membrane thermal fluctuations to this process, via the attractive 
interactions that they mediate between the adhesion points. The origin of the fluctuation-induced 
(Casimir-like) interactions are the restrictions imposed on the membrane thermal fluctuations by 
the adhesion points, and the associated free energy cost which is minimized when the adhesion 
points localize in a cluster. We investigated both the two- and many-body fluctuation-induced 
interactions. For the two-body problem, our analysis reveals that the fluctuation induced pair 
potential is infinitely long-range with a logarithmic dependence on the pair distance. If, in 
addition to the excluded volume interactions, the membrane and the surface also interact via an 
attractive confining potential, the fluctuation-induced pair potential becomes screened at large 
distances. The screening of the pair potential is due to the fact that far away from each adhesion 
point, the amplitude of the fluctuations is governed by the strength of the external potential 
rather than by the presence of the other adhesion point. 

In the many-body problem, the fluctuation-induced interactions are self-screened. The am- 
plitude of the thermal fluctuations at each unit area of the membrane is governed by the distance 
to the closest adhesion points, which implies that each point interacts only with a few nearby 
points. This justifies our mapping of the problem into the 2D lattice-gas model with an effec- 
tively larger ( "renormalized" ) interaction energy. Depending on the strength of the renormalized 
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interactions, the system may be either in a "gas" (uniform distribution) or a "condensed" (ad- 
hesion cluster) phase. The interesting question which arises is whether the fluctuation-induced 
contribution to the attraction is sufficiently strong to allow cluster formation. Our analysis finds 
that the answer to this question is no. The fluctuation-induced interactions alone are too weak 
to induce the condensation transition. They do, however, greatly reduce (to below the thermal 
energy ksT) the strength of the direct interactions at which the transition takes place. 
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